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In this paper we compare the dimensional method with the Lie groups tactic in order to show the 
limitations and advantages of each technique. For this purpose we study in detail a perfect fluid 
cosmological model with time-varying "constants" by using dimensional analysis and the symmetry 
method. We revise our previous conclusion about the variation of the fine structure constant finding 
for example that in the radiation predominance era if a varies is only due to the variation of e 2 e^" 1 
since ch = const in this era. 



PACS numbers: 98.80.Hw, 04.20.Jb, 02.20.Hj, 06.20.Jr 



I. INTRODUCTION 

The main goal of this paper is to show how fruitful is 
to work with dimensional techniques (DA) as compared 
to other more sophisticated methods as the Lie groups 
(LG). For this purpose we study some cosmological mod- 
els that consider all the physical "constants" i.e. G, c, A 
and H as functions depending on time t. In particular we 
will study models which matter content is modeled by 
a perfect fluid. Therefore this paper is devoted to com- 
pare the dimensional method with the Lie group tactic, 
by showing the limitations and advantages of each tech- 
nique. 

We understand that the dimensional method is in ac- 
tually a method of Lie groups, this method has such al- 
gebraic structure, but the Pi-Theorem only finds scaling 
symmetries while the Lie method finds all the possible 
symmetries. One of our purposes in this paper will be to 
show this fact. 

Cosmological models with time-varying "constants" 
have been studied for quite some time ever since Dirac 
proposed a theory with a time- varying gravitational con- 
stant G. Several works have investigated cosmologi- 
cal models with variable cosmological constant within a 
framework of dissipativc thermodynamics as well as in 
the case of perfect fluids. 

The purpose of this work is to perform a detailed study 
of all the possible symmetries of a perfect fluid model 
with time varying constants showing that in this case it 
is possible to find more solutions in addition to the scal- 
ing one (obtained through DA). In order to carry out this 
study, we begin in section 2 by outlining the equations 
that govern the model as well as the notation employed. 
We present three models. The first of them will be for- 
mulated without the condition div(T*) = 0. In order to 
make that "constant" h appears into the field equations 
we impose an adequate equation of state for the energy 
density, in this way through this condition we will be 
able to formulate a very general equation that contains 



to all the "constants" . The second and third of the mod- 
els verify the condition diviTj) = 0. As in the first of our 
models we will consider two cases in each model, to take 
a general form for the energy density and the case which 
verifies the equation of state p = aB 4 (black body radia- 
tion) . In order to get rid of the entropy problem in our 
third model we will consider adiabatic matter creation, 
showing that for an adequate valour of a /3-paramctcr this 
model could be reduced to our second model. Therefore 
this will be our more general model. 

In section 3, we review the scaling solution obtained in 
previous works by highlighting the "assumed" hypothe- 
ses that we need to make in order to obtain a solution 
using dimensional analysis, these are: div(T*) = 0, con- 
servation principle, and that the relation G/c 2 remain 
constant for all values of t (cosmic time). In this section 
we will find a solution for the field equations through two 
different dimensional ways. In the first of them we apply 
the Pi-theorem in order to solve our most general model, 
the third one, and as we will able to see this is, in our 
opinion, the best method to solve the equations since it is 
the simplest one and it lets us to obtain a complete solu- 
tion. We discuss some interesting relationships that arise 
with the similarity method. We will revise our conclu- 
sions obtained in previous works trying to improve them, 
in this way we arrive to the conclusion that the fine struc- 
ture constant varies but we cannot know which constant 
or constants are the responsible of such variation since we 
can only calculate the behaviour of the constants G, c, A 
and h but no of any electromagnetic constant, e,SQ,.... 
With regard to the second dimensional tactic we will try 
to solve the first of our model, which does not verify the 
condition div(TJ-) = 0. In order to do it we shall need 
to impose some hypotheses to reduce the number of un- 
known quatities in the field equations. We arrive to the 
conclusion that the solution obtained verifies the general 
conservation principle (this is a limitation of the dimen- 
sional tactic) and that it is less general than the solution 
obtained with our first dimensional method. This solu- 
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tion is only valid for the case of radiation predominance 
i.e. uj — 1/3 while the above solution is valid for all value 
of uj that is to say for all kind of matter and not only for 
radiation. 

In Section 4, we work towards finding other possible so- 
lutions to the field equations using the Lie group method. 
We start this section by rewriting the field equations in 
such a way that we can use the standard Lie procedure 
that allow us to find more symmetries. With this tactic 
we only will be able the study the second and the third 
model but without the possibility of considering the case 
where p = a9 4 that is to say we are not able to study 
the general case where h appears. After outlining the 
equation and the constraint, we proceed to study some 
cases. The first one is the obtained previously by using 
dimensional analysis since dimensional analysis is just a 
special class of symmetry (scaling symmetry). We would 
like to emphasize that the Lie method show us that one 
of the assumptions made with the dimensional method, 
G/c 2 = const., is at least correct from the mathematical 
point of view. This result allow us to validate completely 
the solution obtained through similarity. Nevertheless 
there are more symmetries, some of them correspond to 
our second and third solutions 

We conclude the paper by summarizing the results in 
section 5 and emphasizing why the dimensional method 
works so well. In the appendix we try to show how te- 
dious is to work with the Lie method if one tries to be 
rigorous. 



II. THE MODEL 

In this section we will outline the field equations of 
three models. In the first of them, its energy-momentum 
tensor does not verify the conservation principle. The 
second model and the third do it but in the third one we 
will consider adiabatic matter creation. In all these cases 
we will consider two possibilities, a general form for the 
energy density p and a specific equation of state for this 
p = ad\ in order to introduce into the field equations the 
Planck constants through the radiation "constant" a. 

We will use in all the models the field equations in the 
form: 



with 



1 8nG(t) 

Rij - ^9tjR = -^TjT-Tij + A(%y, (1) 



where the energy momentum tensor is: 
Tij = (p + p) UiUj - pgij, 



(2) 



and p = ojp in such a way that uj £ ( — 1,1], that is to 
say, our universe is modeled by a perfect fluid. 
The line element is defined by: 



ds 2 



c 2 dt 2 + dn 2 



(3) 



dn 2 = f{t) 



dr 2 



1 - kr 



r 2 (d9 2 + sin 2 



(4) 



we consider "only" a flat model i.e. k — 0. as the most 
recent observations suggest us 1-4 : 
The cosmological equations are now: 

2H' + :iH 2 = -^}^p + c(t) 2 A(t), (5) 



c(t) 2 



3H 2 



inG(t) 



p + c(t) 2 A(t), 



(6) 



where H = (/'//) is the Hubble function. 

Applying the covariance divergence to the second 
member of equation (1) we get: 



(7) 



that simplified is: 



T J ■ 



which yields: 



c G 



8irG 



(8) 



p' + 3(uj + l)pH = - 



a'c 4 a A d 

<37rG G c 



(9) 



where it is noted that this is the main difference with 
respect to other approaches. 

Therefore our first model is governed by the following 
equations: 



2H' + 3H 2 



^p + c(t) 2 A(t), (10) 

c(t) 2 

m 2 = ^M^ P + c(t) 2 A(t), (11) 



P ' + 3(uj + l)pH 



c(t) 2 

A'c 4 G' c! 
~^G- p G +ip -' 



(12) 



In order to incorporate the "Planck constant H" we 
chose or impose an adequate equation of state for the 
energy density, we use the black body equation of state 



7T k 

p = a9 4 where a = 15e3 | 3 so that equation (9) is now 

15A'c 7 ^ 3 



4- -3 



d ti 
~ + J 



+ 3(uj + 1)H = 



S^Gk%9 4 ' G 



c 



(13) 

in this way our first modified model is governed by 
the following equations: 
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2H' + 3H 2 



c(t) 2 



p + c(t) 2 A(t), 



4— 



c 
c 



7T 



3(v + l)H 



15A'c 7 ft 3 G[_ _ d_ 
8-K 3 Gk%9 i + ~G~~c~' 



(14) 
(15) 
(16) 



The second class of models that we study verifies the 
principle of conservation for its energy-momentum tensor 
i.e. we assume that div(Tj) = 0, then equation (9) is 
reduced to: 



P 

'c 4 



A'- 4 



3(oj + l)pH 



G A c 



0. 
0. 



or if p = a8 4 , equation (13) reads now: 

-3(w + l)H 



4- 



87r 3 G7c 4 6> 4 



(17) 
(18) 

(19) 
(20) 



Hence the field equations for our second model are 
(for a general form of p): 



2H' + 3H 2 
3H 2 

p 1 + 3{lu + \) P H 
4< 



8ttG 



Ac' 



c 
8itG 
—P 

0, 
0. 



Ac' 



8npG G 
and for the black body equation of state p 

8itG 



(21) 

(22) 
(23) 

(24) 



a9 4 



2H' + 3H Z 



3H l 
d_ K_ 

c h 

l5A'c 7 h 3 



inG 
-7TP 



-p - Ac 2 = 



Ac 2 



8TT 3 Gk%e i 



3(u> + l)H 

G' c' 
+ G- A ~ 



0, 


(25) 


0, 


(26) 


0. 


(27) 


0. 


(28) 



To end we study briefly the important case in which 
adiabatic matter creation _ 8 can be taken into account, 
in order to get rid of the entropy problem. This will be 
our third model. The matter creation theory is based on 
an interpretation of the matter energy-stress tensor in 
open thermodynamic systems, which leads to the modi- 
fication of the adiabatic energy conservation law and as 
a result including the irreversible matter creation. The 



matter creation corresponds to an irreversible energy flow 
from the gravitational field to the constituents of the par- 
ticles created and this involves the addition of a creation 
pressure p c in the matter energy- momentum tensor which 
we discuss below. 

The field equations that now govern our model are as 
follows: 



2H' + 3H 2 = ~^S^(p + p c ) + c 2 (t)A(t) 7 



'it) 



3H 2 



8nG(t) 



p + c 2 (t)A(t), 



n + 3nH = ip, 



(29) 

(30) 
(31) 



and taking again into account our general assumption on 
the conservation principle i.e. equation (8) with T?j = 0, 
we obtain the two equations 



p' + 3(p + p + Pc )H = 0, 



^T^7^~ 4_ =0 ' 
8nGp G c 



(32) 
(33) 



and where n is the particle number density, ip is the func- 
tion that measures the matter creation, H = /'/ / rep- 
resents the Hubble parameter (/ is the scale factor that 
appears in the metric), p is the thermostatic pressure, p 
is energy density and p c is the pressure that generates 
the matter creation. 

The creation pressure p c depends on the function 
For adiabatic matter creation this pressure takes the fol- 
lowing form: 



Pc = 



P 



3nH 



(34) 



The state equation that we next use is the well-known 
expression p = Lop, where u> = const, and lo G (— 1, 1]. 
Therefore, the new set of field equations is now: 



2H' + 3H 2 = - 



8nG(t) 



3H 



c 2 (t) 
8nG(t) 



(p + Pc ) + c 2 (t)A{t), 



c 2 (t) 

p' + 3(cj + l)pH = (u> + l)p 
A'c 4 



(35) 

P + c 2 {t)A{t), (36) 



9l 

8-kGp ' G c 



= 0. 



(37) 
(38) 
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If we study from the dimensional point of view the 
equations (35-38) it is found the following relationship 
between the quantities: 



7Tl 



7T3 



7T5 



Gpt 2 



Gpt 2 
c 2 
n 



TT 2 



7T 4 



7T6 



A'c 4 



7T 7 



c 2 Ai 2 , 

Gpc' 
: A'c 5 ' 



(39) 



where 7Ti and 7T2 have been obtained from equation (35) 
while 7T3 and 7T4 have been obtained from equations (36), 
and 7T5 from (37). The 7r — monomias ttq and 7T7 have 
been obtained from equation (38) 

In this way we have the following conclusions: 

1. From the monomias w± — 1^3, we see that the quan- 
tities p,p c and p behave in a similar way i.e. 



p»p c & p. 



(40) 



as we already know, since we have imposed the 
equation of state p « p, but the DA of the field 
equations suggests us that p c « p, i.e. that the 
pressure p c behaves as the energy density. 

2. From 114 = c 2 At 2 we obtain a clear behaviour for 
the cosmological constant: 



A 



1 

'cH 2 



3. From tt^ it is obtained that 
tp « nH, 



(41) 



(42) 



i.e., the DA suggest us that ip must be propor- 
tional to the Hubble parameter. Of course other 
possibilities are allowed, see for example the model 
presented by Prigogine and collaborators, where 
ip « H 2 . 



4. And from the monomias ttq and it-? we see that 

9L ~ - 

~G~~' 



(43) 



It is noted that the DA does not understand of numer- 
ical factors, only relations between quantities. 

Since the dimensional method has suggested 
(ip « nH) , we assume that the matter creation function 
follows the law (the same than Lima et al 8 : 



tp = 3(3nH, 



(44) 



where /3 € [0, 1] is a dimcnsionlcss constant (if j3 = then 
there is no matter creation since ip = 0). The generalized 
principle of conservation = 0, for the stress-energy 
tensor (32) leads us to: 



In this way, the set of field equations for the third 
model are now: 



2H' + 3H 2 



(p + Pc )~c 2 (t)A(t)=0, (46) 



8irG{t) 

3H 2 -^p-c 2 (t)A(t) 
pf + 3{u + 1) (1 - 0) pH 

1 4— 

G c 



SnGp 

or equivalently for p = a6 4 : 
8irG(t) 



2H' + 3H 2 



c 2 (t) 
3H 2 - 

h' 



(p + Pc )-c 2 (t)A(t) 



~dW 



\t)A(t) 



+ 3(w + l) (!-#)# 



87r 3 Gfct6» 4 + ~G ~ A ~c~ 



0. 
0. 
0. 

0. 

0. 

0, 
0. 



(47) 
(48) 
(49) 

(50) 
(51) 
(52) 
(53) 



pf + 3(w + 1) (1 - 0) P H = 0. 



(45) 



We must emphasize that these models with (3 = (no 
matter creation) reduces to the second class of models. 

Therefore these are the more general field equations. 
In the next sections we will find a solution for all these 
field equations through different ways beginning with the 
dimensional one and ending with the Lie method. 



III. DIMENSIONAL METHODS. 

In this section we will find a solution to the field equa- 
tions through two different dimensional ways . We 
will begin in the next subsection, the simplest method, 
finding a solution to eqs. (50-53) applying the Pi- 
theorem, while in the following subsection, the not so 
simple method, we will find a solutions to the field eqs. 
(14-16) applying a dimensional tactic which consists in 
reducing the number of variables into the equations mak- 
ing a simple hypotheses obtained through DA. In this way 
it is obtained a simple differential equations that admits 
a trivial integration. 



A. The simplest method. 

Our purpose in this subsection is to integrate eqs. (50- 
53) that is to say to find a solution to these equations. 
In this case we apply the Pi-theorem. In order to apply 
the Pi-theorem in the first place we need to fix the set 
of governing parameters p = p(C,9Jl) which it is com- 
posed of the set of fundamental constants C and the set 
of fundamental quantities 3JI. The set of fundamental 
constants C contains the physical constants and the char- 
acteristic constants of the model. This model evidently 
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has no physical "constants" since these vary. Our first 
characteristic constant will be obtaining by integrating 
equation (45), that brings us to obtain the following re- 
lation between the energy density and the scale factor 
and which is more important, the constant of integration 
that we shall need for our subsequent calculations: 



-3(w+l)(l-/3) 



(54) 



where A Ut p is the integration constant that depends on 
the equation of state that we need to consider i.e. con- 
stant u> and constant (3 that controls the matter creation, 
[A^.p] = L^+W-fi^MT- 2 , where we arc using a 
dimensional base 03 = {L, M, T} , see 15 for details. To 
apply the Pi-theorem we need another constant. This 
new constants is obtained by making a simple hypothe- 
sis about the behaviour of the "constants" G and c. We 
suppose that the relation G/c 2 = B, [B] = LM^T , 
remain constant in spite of both "constants" vary, but 
in such a way that this relationship remain constant for 
all t, the universal time. Furthermore, in this way we 
are guaranteeing that it is verified the covariance prin- 
ciple. With the Lie group tactic we will show that at 
least this relationship has mathematical meaning since 
it is deduced as solution of the field equations and not 
imposed as hypothesis. Therefore the set of fundamental 
constants is: C = C (Aj„3, B) . Now, the set of fundamen- 
tal quantities has only one quantity, the universal time, 
t, as it can be trivially deduced through the Killing equa- 
tions. Therefore 971 = 9JT(t) . 

Our purpose is to show that no more hypothesis are 
necessary to solve the differential equations that govern 
the model. Therefore the set of governing parameters 
is now: p = p{A^^,B,t) , that brings us to obtain the 
next relations: 





G 


Au.p 


B 


t 


L 


3 


1 


1 





M 


-1 


1 


-1 





T 


-2 


-1 





1 



. 2 2 , i 2{1 — -y) 

Goc A^+p B—i +1 t^+^ 



(55) 



In this way, it can be easily obtained the rest of quan- 
tities, obtaining: 

— — i (i — f) 
c oc A^gB-y+H T+ 1 , 

ha: A^ B^tt- 1 ^ , 



miOtA^B 7+if-y+i, 
potB-H- 2 , 



(56) 



3 _1 4 ~ 2 7 



k B 6<* A^B^'h^t, 
a - 1 / i s oc A^B^'h^-%, 
A cx A^B^t^, 

7-1 

1=2 

where 7 = 3(o; + 1)(1 — 0) — 1, and q is the decelera- 
tion parameter. If lu = 1/3 then it is observed that the 



entropy is not constant 

a~ 1/4 s oc t^^' 



(57) 



we can check that the next results are verified: we see 
that § = B\f(3 (trivially), p = aB^ V/3, / = ct V/3, A oc 
pp- oc f~ 2 while the relation he 7^ const, since it depends 
on (3 (if = 0, then He = const but only when oj = 1/3 
i.e. in the radiation predominance era). Wc emphasize 
this relationship between these two "constants" because 
as it is known from the definition of the fine structure 
constants a, 



(58) 



ineofic 

if a varies (in this epoch and only in this one) this vari- 
ation only can be caused by e 2 /eq- While in other epoch 
he =/= const independently of (3, being very difficult to ex- 
plain the origin of such variation or more exactly which 
one or which ones of the constant are the cause of such 
variation. With this model we can only calculate the vari- 
ation of the "constants" c and h (and if we are strictly 
rigorous we would have to say that only for the time of 
radiation predominance) and we cannot calculate the be- 
haviour of any electromagnetic quantity. Wc believe that 
this must be the behaviour between this "constants" be- 
cause as we have shown in a previous paper"' other rela- 
tion as h « c brings us to a static universe. Furthermore, 
the behavior obtained here works well in the framework 
of the quantum cosmology as we have pointed out in 
reference 17 . 

In previous works we calculate the behaviour of e 2 E§ 1 
as a function of (Aj,/3, B, t) in such a way that a (dimcn- 
sionlcss quantity) remain constant in spite of the fact 
of that all "constants" vary but in a conspire way since 
they vary but keep alpha constant. Now we believe that 
the relationship, for e 2 e^ 1 was wrong in spite of the fact 
that there is a relation between c and Eq, c 2 — (/io^o) 1 
. We cannot obtain an expression for electromagnetic 
quantities only in function of (A u ^,B,t) (obtained in 
the framework of standard cosmology). This fact has 
been pointed out in a previous work (but in a differ- 
ent framework since in that work c = const.) where we 
showed that it is impossible to reconcile the cosmologi- 
cal quantities with the electromagnetic quantities (see 18 
for details). Here we have the same situation. Therefore 
alpha can vary but we do not know which of the con- 
stants are the responsible of such variation as already it 
has been pointed out by M. E. Tobar 19 . 

We also can check that our model has no the so called 
Planck's problem since the Planck system behaves now 
as: 



m 
t 



p Kg) 



fit), 



(59) 



since the radius of the Universe f(t) at Planck's epoch 
coincides with the Planck's length f(t p ) « l p , while 
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the energy density at Planck's epoch coincides with the 
Planck's energy density p(t p ) m p p w t~ 2 , where p p = 
m p c 2 /l p . See 1 ' for more details and the followed method 
etc... 

It is observed from (56) that if we make (3 = the 
following set of solutions are obtained: 



UJ 


1 


1/6 


1 /Q 

1/3 


U 


1 /Q 

-1/3 


— Z/6 


f 


1/3 


2/5 


1/2 


2/3 


1 


2 


P 


-2 


-2 


-2 


-2 


-2 


-2 


6 


-1 


-4/5 


-1/2 





1 


4 


s 


-1/2 


-3/10 





1/2 


3/2 


9/2 


G 


-4/3 


-6/5 


-1 


-2/3 





2 


c 


-2/3 


-3/5 


-1/2 


-1/3 





1 


h 





1/5 


1/2 


1 


2 


5 


rrii 


1/3 


2/5 


1/2 


2/3 


1 


2 


A 


-2/3 


-4/5 


-1 


-4/3 


-2 


-4 


q 


5/2 


3/2 


1 


1/2 





-1/2 



with: (uj = — 1) corresponds to de Sitter (false vacuum) 
represented by the cosmological constant (special case), 
(uj = — |) for domain walls, (u = — -|) for strings, (u> = 0) 
for dust (matter predominance), (uj — ^) for radiation 
or ultrarelativistic gases (radiation predominance), (uj = 
|) for perfect gases, (uj = 1) for ultra-stiff matter. For 
example, if we take the case uj = the table (60) tells us 
that 

/ oc i 2/3 , pcx£~ 2 , oct° = const., 

G cx i _2/3 , coci~ 1/3 , h oc t, ...etc. (61) 

This table tells us too that if we want that our universe 
accelerates then we have to impose that — 1 < ui < — 1/3 
but we must be careful since with this parameter we see 
that the temperature increases. 

We can try to generalize this scenario taking into ac- 
count various kinds of matter. The idea is as follow. We 
can define a general energy density p as: 

6 

P = X)ft (62) 

where pi stands for each kind of energy density, and the 
parameter i = 0, 1, 6 in such a way that: i = corre- 
spond to uj = —1 (the false vacuum), i = 1 correspond 
to domain walls i.e. to u> = —2/3, i = 2 to w = —1/3, 
i = 3 to uj = 0, i = 4 to uj = 1/3, i = 5 to uj = 2/3 and 
finally i = 6 to u> = 1; but in such a way that each type 
of matter verifies the relation 

Pi = ^tPi (63) 
in this way we define the total pressure as: 

6 

p = Y,Pi (64) 
«=o 

but in this case we do not impose that it is verified for 
each kind of matter the relation: 

Pi =A u J- 3 ^ (65) 



Our purpose is as follows: we impose that the relation 
div(T) = it is verified i.e. 

p' + 3(p + p)H = (66) 

taken into account that pi = uiipi then 

6 

f? + 3Hj2[(uJi + l)pi]=0 (67) 

that we can "approximate" through the relation: 

p = A m .r m (68) 

with m = 3J^ =0 (uji + 1) . It is proven in a trivial way 
that if we consider only one type of matter we then re- 
cuperate the above results i.e.. m = 3(cj + 1). 

Therefore with the next set of governing quantities 
p = p (A m , B, t) we arrive to obtain the following table 
of results: 

G oc A*^ B^t^~ 2 
cocA^B^t^' 1 

-2 _ 2 _ 4 

A (X Am 1 B^+it^+1 

hocA^B^TT-H^TT- 1 

-J— -x 2 

rrii oc Am 1 B*+ 1 t :c + 1 (69) 

f (X Am 1 B~t~ 

pot B~H- 2 

k B e ex a^b^tt- 1 ^- 2 

a" 1 / 4 s cx A^B^tt-^^tt-I 

where x + 1 = m . = 3 X)i=o ( Wi + 1) ■ ^ f° r example 
we consider an universe with dust (uj = 0) and radiation 
(uj = 1/3) then m = 9, wc obtain a very surprising results 
as we can see 

G « T 8 / 5 , ccx^ 4 / 5 , ficxt" 2 / 5 , Acxi- 2 / 5 , 
p (x i~ 2 , /cxi 1/5 etc... (70) 

As we have seen, it is obtained a complete solution 
of the field equations making only one hypothesis i.e. 
G/c 2 = B. The method, a direct application of the Pi- 
theorem, allows us to obtain, in a trivial way, the be- 
haviour of all the quantities under study. Furthermore, 
the followed tactic allows us to generalize the model con- 
sidering a mixed of fluids (but without interactions be- 
tween them). This method is therefore, simple but pow- 
erful as we will see in the next sections since with the rest 
of the tactics we will not be able to solve this so complex 
model. 

B. Not so simple method. 

In this subsection we will try to solve the first of our 
models eqs. (14-16), i.e. a model where all its "con- 
stants" are time- varying and its energy-momentum ten- 
sor does not verify the conservation principle div(Tf) ^ 
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0. In particular we are interested in solving the eq. (16) 
which contains all the information about the model. As 
this equation has 6 unknown quatities, we will need to 
make some hypotheses in order to simplify the original 
equation and to try to integrate it. In this case we will im- 
pose some hypotheses about the behavior between some 
of the quantities through DA in such a way that these 
hypotheses allows us to reduce the number of variables 
in the equation under study. 

To solve equation (16) we imposed two simplifying 



hypotheses, the first one, that the relation G/c 2 = B 
(where B is a const.) remains constant, and the sec- 
ond one, that the cosmological "constant" verifies the 
relation A oc with d £ K (this is a very strong 

condition), while div(T?) ^ 0, in this way the equation 
was solved perfectly. Since the constant B has dimen- 
sions [B] = LAI^ 1 we can get the dimensionless mono- 
mia 7Ti = PSr- where b £ R. With these hypotheses and 
tti equation (16) simplifies to: 



J 



4- - 3 



T 



■3(u + l)H- 



15d c 4 [c't + c] h 3 

4^3 Gkien 3 



(71) 



that has no immediate integration. We have to take into 
account the field equations (15) 



3H 2 = 8wbr 



dt~ 



(72) 



from this we get / = KJ,* where x = (»2±W) 2 
substituting in (71) together with G = Be 2 we get 



and 



4- - 3 



3(u + l)x 15d c 2 [c't + c] h 3 d _ 
t 4^3 Bk%6H 3 2 ~~ 0, 



r 



(73) 



i.e. one equation with 3 unknowns. 

If we want to integrate equation (73) we have to take 
a decision on the behavior of the constant h. 



Taking ch = const. « K = — then j- = — ^ yielding 



4- 



3(w + l)x Ihd A 3 [c't + c] 



4tt 3 cBk%6H 3 



2°- 

c 



(74) 



Also if p = a0 A and p = 



b 

~W 2 



b 

"BT 3 



a6 we have: 



f i5c 3 (t)n 3 (t)b y (i5A 3 by 1/2 



(75) 



And substituting into the previous equation we get: 
-2 3{cj + 1)k 15dA 3 [c't + c] 



4tt 3 cBt 3 ^^ 



2°- 

c 



0. 



t t 

and simplifying 

-2 3{u + l)x d 
T H t 4?r6 

therefore we get a very simple differential equation: 



(76) 



d 1 
- + 1 



2 C - = 0, 

c 



(77) 



c 
c 



127r&(w + l)^-87rfr-cZ 



8nb + d 



(78) 



integrating it we obtain easily: 
c = K s tt 



127rb(^ + l)>t-87rb-rf 
8-nb+d 



(79) 



where £ 

We can consider another possibility. Take the group 
of governing quantities p=p{K^, A,t} where is the 
proportionality constant obtained from / = K^t* and A 
is the constant establishing the relation between h and 
c. The results obtained by means of the gauge relations 
are: 

G oc K^A-H 6 *- 4 
cozK^t*- 1 
hoc K- 1 At 1 -* 

k B e oc K-^ASr" (80) 
p oc K- 4 At- 4 * 
rrii oc KlAt.- 3>t + 2 
A oc K-H- 2 * 

where mj comes from the energy density definition pe = 

2 

""i c (n stands for the particles number) . We can check 
that we recover the general covariance property 5? = 
const, iif yc = \. Similarly we can see that the following 
relations are satisfied: p = a6 4 , p = Af~ 4 (equivalent to 
div(T?) = 0), A oc f~ 2 and f = ct (no horizon problem). 
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With the value h — \ we get 



coci- 1 / 2 , next 1 / 2 , Gcxt^ 1 , k B 0(xt-^ 2 , 01 , 

-2 ™.^ t l/2 ( 81 ) 



To obtain this solution we have needed three hypothe- 
ses: 

1. G/c 2 = B, 

2. A w i~ 2 , 

3. cfi, = const. 

Some of them are very restrictive since they induce 
a scaling solution (power law solution) and as we can 
see we have obtained the same solution than in the pre- 
vious subsection for a model with w = 1/3 and (3 = 0, 
i.e. with radiation predominance and without matter cre- 
ation. But. How it is possible?. We have tried to solve 
eq. (16) i.e. a model that, in principle, it does not verify 
the condition divt{Tf) = however our solution verifies 
such condition. As we will see in the last section, the 
DA always is related to conservation principles, for this 
reason in spite of working with eq. (16) our hypotheses 
obtained through DA have taken us to solve our second 
model which verifies the condition div(T?) = 0. 

Furthermore, we must emphasize that with this tactic 
we have obtained less information since our solution is 
only valid for the case to = 1/3, while with the naive 
method, Pi-theorem, the solution is valid for all kind of 
matter. 

Therefore our first dimensional method, the simplest 
one, needs fewer simplifying hypotheses and its solution 
is much more general that the obtained one with this 
other dimensional method. 



most general form of symmetries. Hence, we shall study 
the model through the method of Lie group symmetries, 
showing that under the assumed hypotheses there are 
other solutions of the field equations. In this section we 
shall show how the lie method allows us to obtain dif- 
ferent solutions for the field equations. In particular we 
seek the forms of G and c for which our field equations 
admit symmetries i.e. are intcgrable (see 21-30 ). 

An alternative use of the Lie groups have been per- 
formed by M. Szydlowki et. al. 51-32 where they study 
the Friedman equations in order to find the correct equa- 
tion of state following pioneer works of Collins 33 . 

In order to use the Lie method, we rewrite the field 
equations as follows. From (21) — (22), we obtain 



and therefore 



2(fl)' = -^(p + p). 



From equation (23), we can obtain 



H = 



3((w + l)p)' 



therefore 



= 12vr(u; + l) 2 ^L. 
P J c 



(82) 



(83) 



(84) 



(85) 



Taking 12w (cj + 1) = A and then expanding, we obtain 



IV. LIE METHOD 

As we have seen earlier, the ir — monomia is the main 
object in dimensional analysis. It may be defined as a 
product of quantities which are invariant under changes 
of fundamental units, ir — monomia are dimensionlcss 
quantities, their dimensions are equal to unity. Dimen- 
sional analysis has the structure of a Lie group 2 ". The 
7r — monomia are invariant under the action of the simi- 
larity group. On the other hand, we must mention that 
the similarity group is only a special class of the mother 
group of all symmetries that can be obtained using the 
Lie method. For this reason, when one uses dimensional 
analysis, only one of the possible solutions to the problem 
is obtained. 

As we have been able to find a solution through dimen- 
sional analysis, it is possible that there are other symme- 
tries of the model, since dimensional analysis is a reminis- 
cent of scaling symmetries, which obviously are not the 



p" = P —+A%p\ (86) 

p cr 

If we consider the case in which there is matter cre- 
ation, the resulting equation to study is now: 

P "= P -+A%p\ (87) 
p &■ 

with 12tt (lu + if (1 - (3) 2 = A. Therefore, we obtain the 
same equation that we have obtained, namely equation 
(86) except the constant A, that incorporates all the pa- 
rameters that controls the matter creation. 

Now, we apply the standard Lie procedure to this equa- 
tion. A vector field X 

X = t(t,p)d t +r)(t,p)d p , (88) 

is a symmetry of (86) iff 
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<ft - vf P + Vtt + (2vtp - ttt) p' + (Vpp - 2£t P ) P' 2 - ZppP' 3 + 



+ ( ??P - 2£ t - ip'i P ) f - [ V t + (v P - ft) ft - p' 2 Q U = °- 

I 



(89) 



By expanding and separating (89) with respect to powers of p 1 , we obtain the overdetermined system: 
I 



r\tt 



A 



q_ 

,■2 



2G- 



V P p-2£tp + P V-P Vp = 
G 

2r)t P ~ 6t - 3A—p 2 £, p - 2p- 1 i lt = 



G 



G 



p't - 2 V A- P + fa, - 26) A-p 



0. 
0. 

0, 
0. 



(90) 
(91) 

(92) 
(93) 



r 



Solving (90-93), we find that 

£(t,p) = -2et + a, r]{t, p) = (bt + d) p, (94) 
subject to the constrain 

bt + d - 4e 



^=2 C - 
G c 



2et-a 



(95) 



with a, b, e, and d as constants. In order to solve (95), we 
consider the following cases. 



A. Case I: b = and d - 4e = 



In this case, the solution (95) reduces to 



G" 



G 



— = 2 — — = B = const. 
G c c 2 



(96) 



which means that "constants" G and c vary but in such 
a way that the relation ^ remains constant. 

The solution obtained through Dimensional Analysis 
needs to make this relations as hypothesis in order to 
obtain a complete solution for the field equations. This 
case shows us that such hypothesis is correct (at least 
has mathematical sense). 

The knowledge of one symmetry X might suggest the 
form of a particular solution as an invariant of the oper- 
ator X i.e. the solution of 



dt 



dp 



(97) 



this particular solution is known as an invariant solution 
(generalization of similarity solution), therefore the en- 



ergy density is obtained as 

dt dp 
—2et + a Aep 

for simplicity we adopt 



1 



(2et-aY 



P = Pat , 



(98) 



(99) 



Once we have obtained p, we can obtain / (the scale 
factor) from 



Au f-^+l) => f= (^t)**!! 



(100) 



in this way we find H and from eq. (22), we obtain the 
behaviour of A as: 



c 2 A = 3H 2 - 



8irG 

— o-P> 



and therefore, 



A=(3f?-8nB P0 ) ~^ = ~^- 



(101) 



(102) 



If we replace all these results into eq. (24), then we obtain 
the exact behaviour for c, i.e., 



1 c'\ c' 
- + - A= - 

t c I c 



(103) 



where A 



87T g po , with A € R + , i.e. is a positive real 
number and thus, 



c = c t a , 



(104) 



( TTa) 
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Hence, in this case we have found that (see fig.l): 

G = G t- 2a ,c = c t- a ,A = A t- 2( - 1 - a \ 

f = (A u t)^ ,p = Po t~ 2 . (105) 

This is the solution that we have obtained with dimen- 
sional analysis in the previous section. 

B. Case II, b = a = 



In this case, we find that 
G 



Bi> 



(106) 



where h = S — 2 and 6 = On following the same 
procedure as above, we find that 



dt dp s 
-r = — P = Pot , 



(107) 



we must impose the condition sign(d) = sign(e), i.e., 
5 G R + , in order that the solution has some physical 
meaning that the energy density is a decreasing function 
of time t. It is observed that if d — 4e then we obtain 
same solution that the obtained one in the case I. The 
scale factor is found to be 



/ = Kjt 



(108) 



where Kf is an integration constant, and therefore, the 
Hubble parameter is: 



H = 



5 



3(w + i)r 



(109) 



which is similar to the scale factor obtained in case I. To 
obtain the behaviour of the "constants" G, c and A, we 
follow the same steps as in case I, i.e., from 



1 A = 3H 2 - 



8ttG 



we obtain the behaviour of A being: 

/ 



A = 

where, I = (Ki - K 2 ), K x ■■ 
i.e., l6l+. Therefore, 



(110) 



(111) 



3(w+l) 



q-pp- and K 2 = 8irp B 



A' 



21 



(112) 



If we substitute all this results into the next equation 
Ac 4 G' 



inGp G 



0, 



(113) 



we obtain an ODE for c, i.e., 



-(X-2) = -(X-2 + S)- 
c t 



where, A 



47rpo B 



A E 



c = cot 



which leads to 



(114) 



(115) 



with a = ^1 + j^^j such that a e [0, 1). In this way we 
can find the rest of quantities: 



G = G t- 2{a+1)+ " 



A = An<- 2(1 - Q) , 



(116) 



note that a < 1. The case a = 1 5 = is forbiden 
and a = brings us to the limiting case of the G, A 
variable cosmologies 41 . 

We notice that this solution is very similar to the case 
I but in this case all the parameters are perturbed by 5 
and more important is the result, ^ = Bt^ (sec figs. 2, 3 
and 4). 



C. Case III, b = e = 

Following the same procedure as above, we find in this 
case that such restrictions imply £(t, p) = a, r](t, p) = dp 
and therefore: 



- = 2---, 

G c a 



which brings us to: 



G 



K exp(— at), 



(117) 



(118) 



where - = a and note that [A'] = [B] i.e has the same 
dimensional equation, 



dt 



dp 



dt dp 

t(t>P) v(t,P) a dp 



p = Pa exp (at) 



(119) 

this expression only has sense if a G R , note that [a] = 
The scale factor / satisfies the relationship: 

p = A^f- Z ^ +1) =>f = K f exp (crf)^TT> , (120) 

that is to say, it is a growing function without singularity. 
In this way, we find that 



H= — 



= cants. 



H>0. 



3(w + l) 

The cosmological "constant" is obtained as 



(121) 



c 2 A 



3(w + iy 



8-n-Kpa =^ c 2 A = I, (122) 
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FIG. 2: Time variation of c(i) for the second class of solutions 
for different values of a : a = 0.5 (solid curve), a = 0.9 
(dotted curve), a = 0.1 (dashed curve) and a = 0.000001 
(long dashed curve), the last solution describes the case c(t) = 
const. 




FIG. 3: The variation of the gravitational "constant" G(t), 
for different values of a and 5. : a — 0.5 and 5 = 1 (solid 
curve), a = 0.9 and 8 — 1 (dotted curve), a = 0.1 and 5 = 1 
(dashed curve) and a = 0.000001 and 5 — 5 (long dashed 
curve), the last curve describes a growing solution. 



note that [I] = T~ 2 , if we replace all these results into 
eq. (24) then we shall obtain the exact behaviour for c, 
i.e. 

(-^— + 2) C - = a, (123) 
\8nhpo J c 



FIG. 4: Time variation of A(t) for the second class of solutions 
for different values of a : a = 0.5 (solid curve), a = 0.9 
(dotted curve), a = 0.1 (dashed curve) and a = 0.000001 
(long dashed curve). In all cases, A(i) is a decreasing function. 

and hence, 

c=Kexp(cot), (124) 
where c = — ^ r- with Co € M _ since a € that 

{ SnKpo + 2 ) 

is, c is a decreasing function on time t. 
In this case, we have found 

c = ifexp(c i), (125) 
G = G exp((-a + 2c )<), (126) 
A = Zcxp(c ir 2 , (127) 

therefore the solutions for this case are (see fig. 5): 

G = G exp((— a + 2c ) t), c = Kexp(c t), 
A = I exp(cot)~ 2 , p = po exp (at) , 

f = K f cxp(at)*^V . (128) 

We have only obtained three solutions but playing with 
the constants a,b,d and e more solutions can be obtained 
(without forgetting that we are only interested in solu- 
tions with physical sense). As we indicated in the intro- 
duction and throughout all the paper, we have obtained 
in the first of our solutions (Case I), the same solution 
that the obtained one with the dimensional methods. In 
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FIG. 5: We see the behaviour of constants for the third class of solutions for different values of cq : cq = —0.5 (solid curve), 
co = —0.3 (dotted curve), co = —0.01 (dashed curve) and a — —0.1. In all cases, A(t) is a growing function. 



this occasion we have not needed to make a previous hy- 
pothesis in order to arrive a complete solution, this is 
therefore one of the advantages of this method as op- 
posed to the dimensional one. We have also studied two 
other cases which can be considered as physically rele- 
vant solutions since / is a growing function on time and 
p is a decreasing function on time. They could describe 
very early cosmological solutions (inflationary ones) . 

We see that with this sophisticated method we can 
only solve the models 2 and 3 in which the energy den- 
sity is expressed in a generic form i.e. we cannot con- 
sider the complicated case in which appears the Planck's 
constant H (the disadvantage). Although we have not 
been excessively scrupulous (formal in the procedure (see 
appendix)) the method becomes more and more compli- 
cated in function of the complexity of the equation under 
study. 

V. WHY THE DIMENSIONAL METHOD IS SO 
WONDERFUL AND CONCLUSIONS. 

As have been pointed out by Carr and Coley 34 , 
the existence of self- similar solutions (Barenblatt and 
Zcldovich 3 5 ) is related to conservation laws and to the 
invariance of the problem with respect to the group of 
similarity transformations of quantities with independent 
dimensions. This can be characterized within general rel- 
ativity by the existence of a homotetic vector field and 
for this reason one must distinguish between geometrical 
and physical self-similarity. Geometrical similarity is a 
property of the spacetime metric, whereas physical sim- 
ilarity is a property of the matter fields (our case). In 
the case of perfect fluid solutions admitting a homotetic 
vector, geometrical self-similarity implies physical self- 
similarity. 

As we show in this section as well as in previous works, 
the assumption of self-similarity reduces the mathemat- 
ical complexity of the governing differential equations. 
This makes such solutions easier to study mathemat- 



ically. Indeed self-similarity in the broadest Lie sense 
refers to an invariance which allows such a reduction. 

Perfect fluid space-times admitting a homotetic vector 
within general relativity have been studied by Eardley . 
In such space-times, all physical transformations occur 
according to their respective dimensions, in such a way 
that geometric and physical self-similarity coincide. It is 
said that these space-times admit a transitive similarity 
group and space-times admitting a non-trivial similar- 
ity group are called self-similar. Our model i.e. a flat 
FRW model with a perfect fluid stress-energy tensor has 
this property and as already have been pointed out by 
Wainwright , this model has a power law solution. 

Under the action of a similarity group, each physical 
quantity <p transforms according to its dimension q under 
the scale transformation. For space-times with a transi- 
tive similarity group, dimcnsionlcss quantities are there- 
fore spacetime constants. This implies that the ratio of 
the pressure of the energy density is constant so that the 
only possible equation of state is the usual one in cosmol- 
ogy i.e. p = Lop, where ui is a constant. In the same way, 
the existence of homotetic vector implies the existence of 
conserved quantities. 

In this paper we have studied the behaviours of 
time- varying "constants" G, c and A in a perfect fluid 
model. We began reviewing the scaling solution obtained 
through dimensional analysis. 

To obtain this solution, we imposed the assumption, 
div(Tj) = 0, from which we obtained the dimensional 
constant that relates p oc J- 3 ( W +!) an d the relation- 
ship G /c 2 = const. = B remaining constants for all value 
of i, i.e. G and c vary but in such a way that G/c 2 remain 
constant. With these two hypothesis, we have obtained 
a scaling solution for all the quantities. In this context, 
the solution obtained through dimensional analysis show 
us that the "constants" G, c and A are decreasing func- 
tions of time, but in this case decrease slowly than in the 
radiation predominance era, while p and / behave as in 
the FRW model solving the horizon problem. 
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Therefore the DA is very simple and allows us to study 
very complex models. Nevertheless it has very serious 
limitations, we are restricted to models that verify the 
conservation principle div(T?) = and work only with 
flat models i.e. k = 0, for models with k 7^ there are 
not scaling solutions. 

Since we have been able to found a solution through 
similarity, i.e. through dimensional analysis, it is possi- 
ble that there are other symmetries of the model, since 
dimensional analysis is a reminiscent of scaling symme- 
tries, which obviously are not the most general form of 
symmetries. Therefore, we studied the model through 
the method of Lie group symmetries, showing that un- 
der the assumed hypotheses, there are other solutions of 
the field equations (an advantage of this method). 

The first solution obtained is the already obtained one 
through similarity, but, in this case we have showed the 
condition G/c 2 arises as a result and not as an ad- hoc 
condition. We also have studied two other cases which 
can be considered as physically relevant solutions since 
/ is a growing function on time and p is a decreasing 
function on time. They could describe very early cosmo- 
logical solutions (inflationary ones). 

The lie method maybe is the most powerful but has 
drawbacks, it is very complicate. We cannot, we do not 
know, solve the most general of our models, but it allows 
us to find more solutions than with de DA method and 
without the necessity of making previous hypotheses 
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APPENDIX A: COMPLETE LIE METHOD FOR 
THE CASE I (IV A) 

In this section we try to show how the Lie method 
works in order to obtain a complete solution to equation 



n' 2 

p"= P - + Kp\ 
P 



(Al) 



taking K = AB with 12tt (uj + if = A and B = G/c 2 = 
const. We will follow the ways of the canonical coordi- 
nates. 

The first step is to determine the admissible algebra 
L 2 . As in section (IV A) the standard procedure brings 
us to obtain the overdetermined system: 



tpp + p~ x £p 


= 0, 


(A2) 


71 PP -2£_ tp + p- 2 j] - p- l <q p 


= 0, 


(A3) 


2Vt P ~ &t - 3Kp% - 2p- 1 j lt 


= 0, 


(A4) 


rnt-2r]Kp+(ri P -2^)Kp 2 


= 0. 


(A5) 



Solving (A2-A5), we find that 

£(t, p) = b + at, v[t,p) 



-2ap, 



(A6) 



where a is a constant. Thus equation (Al) admits two 
linearly independent operators 



X x = d t , X 2 =td t -2pd p 



(A7) 



The second step consists of determining the type of 
algebra that form these two operators. We see that 

[X lt X 2 ] = X x with £ im - £ 2 7 ?1 = -2p £ 0, (A8) 

therefore X\ and X2 spand a solvable non abelian algebra 
type III (sec Ibragimov (1999) 2 ' for details, theorem 12.6. 
page. 287). 

The next step consists in determining the integrating 
change of variable. Upon introducing the canonical vari- 
ables for Xi (Xis = 0,Xiu = 1) given by 

s = p and u(s) = t, (A9) 

and therefore 

t = u(s) p = s. (A10) 

We transform the operator X\ and X 2 to the form: 

X l = d u , X 2 = -2sd s + ud u . (All) 

their difference from the corresponding operators of type 
III by the factor —2 in X 2 does no hinder the integration. 

We rewrite equation (Al) in the new variables (A9) 
obtaining: 



(1 + Ks 



(A12) 



the canonical variables allow us to reduce the order of 
our ode. If we follow this tactic, then it is obtained the 
next first order ode (type Bernoulli), 



z' = ---Kz 3 s 2 , 
s 



(A13) 



where s = p and z = du/ds = (1/p') and which solution 
is: 



z(s) = ±- 



1 



V2Ks + 6\ ' 



(A14) 



The last step consists in obtaining the solution in the 
original variables, being this: 



p' = p^/2Kp + C 1 
this last ode has the following solution 



2 a ret an 



t 



2p+C 1 
Ci 



c 2 = 



(A15) 



(A16) 
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and therefore 



P 



tan(*±f 



2KC\ 



(A17) 



Therefore this is the most general solution to equation 
(Al) but as it can be observed this seems to be non- 
physical, for this reason we will seek more change of vari- 
ables in order to find a better (physical) solutions. 

We would like to emphasize how the Lie procedure 
brings us to solve equation (A13). Of course this ode 
can be trivially integrated but, at this point we prefer 
to integrate it following two tedious ways. In the first of 
them we study equation (A13) through the Lie method 
i.e. following all the procedure shown above and with the 
second method we try to find and integration factor (Lie 
integration factor) which brings us to a direct integration 
of our equation. 

Therefore studying equation (A13) through the Lie 
method it is obtained the overdetermined system: 



vs + (vz - e.) ( -- - **vj - & (-- - /^v ) 



£(i__2Xz 3 S ) - V (-- s -3Kz 2 sA = (A18) 



Solving (A18), we find that equation (A13) admits the 
following linearly independent operators 



X l = z A s*d z , X 2 = (z - 2Kz i s 6 ) d z , 

z 3 

X 3 = d s d z , X A = sd s --zd z . 

s I 



(A19) 



With one of these fields, we obtain a new set of vari- 
ables which brings us to obtain a new ode which will be 
trivially integrated. For example, from X 3 it is obtained 

i = z(s)s, and h(i) = s, (A20) 
and the inverse transformation is: 



Our second way consists in obtaining a Lie integration 
factor \i (for our ode (A13)) in such a way that p, * ode 
is an exact equation. We find the following integration 
factor: 



1 



which bring us to the following solution 
z = ±- 1 



■y/Ws + Ci ' 



(A25) 



(A26) 



Now we see other changes of variables induced by other 
operators. If for example we calculate the canonical vari- 
ables through the second operator (A7) Xi = tdt — 2pd p 
then 



s = pt 2 and u(s) = hit, 



and therefore 



t = e u{s) and 



1 



P 



(A27) 



(A28) 



e 2u{s) ' 

with these new variables, the original equation yields 

u' (1 + su' + (Ks - 2) sV 2 ) 
u" = 1 — '-, (A29) 

s 

which results more complicate than the original one. Us- 
ing the canonical coordinates it is reduced the order in 
the following way i.e. 



. z(l + zs+ (Ks - 2)sV) 
z = . 



which is first-order ode type Abel, where z = ^ , 

1 



s = pt and z = 



t 2 (p't + 2p) ' 



(A30) 



(A31) 



The solution to eq. (A30) is: 



V-4:Z 2 s 2 + 2Ks 3 z 2 + 4zs - 1 



s = h(i) and z(s) = -r-—, (A21) 

h[i) 



in these new variables equation (A13) yields 

-1 



ti 



Ki 3 ' 



which is integrated by quadratures 



h = 



1 



2Ki 2 



(A22) 



(A23) 



Now we obtain the solution of the equation in the orig- 
inal variables, finding that 



z = ±- 



1 

■y/2Ks + C\ 



(A24) 



+2 arctan 



1 - 2zs 



= 0. 



a/-4z 2 s 2 + 2Ks 3 z 2 + Azs - 1, 

(A32) 

To end we will see another transformation. We are 
interesting into knowing the transformation that pro- 
duces the operator X = tdt + pd p - Knowing that the 
canonical variables are obtained by solving the equations 
(Xs = 0, Xu = 1) we find that 



s = — and u(s) = \nt, 
t w ' 



and therefore 



t = e u{s) and 



p = se 



u(s) 



(A33) 



(A34) 
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with these new variables, equation (Al) simplifying 
yields 



u 
s 

which has a particular solution 

1 1 ( 2 
u = — In 

3 V Ks 



= -u' 2 - - - ( 1 + Ke 3u< - S A su' 3 (A35) 



(A36) 



finding in this way that in the original variables it yields 



2 
Kt 2 



(A37) 



which is the solution obtained in section (IV A) . 

Once again if one insist in solving equation (Al) 
through DA (applying the Pi-theorem, see 38 for details) 
it is founded that with respect to the dimensional base 
23 = {p, T} each quantity has the following dimensional 
equation [p] = p, [t] = t and [K] = p~ x t~ 2 . Therefore, 



we find in a trivial way that: 





P 


K 


t 


p 


1 


-1 





T 





-2 


1 



(A38) 



It is observed that if for example we work with a new 
dimensional base like 23 = {L, M, T} , i.e. the usual one, 
the above quantities have the following new dimensional 
equations: [p] = L- x MT- 2 ,[t] = T and [K] = LM~\ 
Therefore, we have a new overdetermined system of equa- 
tions which determine p, 





P 


K 


t 


L 


-1 


1 





M 


1 


-1 





T 


-2 





1 



As it is observed we have obtained the same solution 
than the obtained one in the previous sections, but now 
making only a simple dimensional considerations. 
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